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ABSTRACT 

To overcome difficulties in understanding the origin of the submillisecond opti- 
cal pulses from SN 1987A we apply a model similar to that of Kundt and Krotscheck 
for pulsed synchrotron emission from the Crab. The interaction of the expected ul- 
trarelativistic e ± pulsar wind with the pulsar dipole electromagnetic wave reflected 
from the w alls of a w pulsar cavity” within the SN 1987A nebula can generate pulsed 
optical emission with efficiency at most T? mwc ~ 10 -3 . The maximum luminosity of 
the source is reproduced and other observational constraints can be satisfied for an 
average wind energy flow « 10 38 erg/(s;steradian) and for electron Lorentz factor 
7 w 10 5 . This model applied to the Crab yields pulsations of much lower luminosity 
and frequency. 
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1. Introduction 


The strong luminosity (between 400 and 900 nm) and the short period (.P — 0.5 ms) 
of the reported optical pulsations from the young supernova remnant (SNR) SN 
1987A (Kristian et al. 1989) raises problems for conventional models of pulsar opti- 
cal emission. If relativistic beaming plays no dominant role, a rather small radiating 
area ^ (ci 3 ) 2 is implied, leading to an extraordinarily high optical brightness tem- 
perature (JfcTt >> 1 GeV). It has not been demonstrated how such emission may 
arise close to a neutron star. On the other hand, it is widely accepted that pul- 
sars may give rise to a wind of relativistic electrons and/or positrons (e*) (Rees 
and Gunn 1974, Kundt and Krotscheck 1977, Kennel and Coroniti 1984, Cheng, 
Ho and Ruderman 1986). As suggested by Kundt and Krotscheck for the Crab 
nebula, ultra relativistic may give rise to pulsed emission far from the stellar 
surface where the relativistic wind runs into the pulsar dipole electromagnetic wave 
reflected from the inner boundary of the surrounding nebula. The main point of 
our paper is that such a mechanism can account successfully for the periodicity of 
the modulated optical signal reported from SN 1987A and it alleviates the optical 
luminosity problem posed by observations. 

During the January 18 observation the brightness of the detected pulsed signal 
varied from magnitude 17 to 16 reaching at its maximum 1% of the luminosity of 
the SN 1987A remnant (Middleditch 1989). Thus, the maximum “optical” pulsed 
luminosity of the source was X 0 pt = 3 • 10 36 erg s -1 x Afl/47r, where A Q, is the 
solid angle into which the pulsed radiation was beamed. At the same time the 
luminosity of the remnant (SNR) was Xsnr = 3 • 10 38 erg/s (Burki and Cramer 
1989). Subsequent observations failed to detect the pulses at a limiting magnitude 
lower by 2 than the maximum observed (Kristian et al. 1989) and by 8 than that of 
the SNR (Ogelman et al. 1989). By the end of April 1989 the remnant bolometric 
luminosity decreased to Xsnr = 1 • 10 38 erg/s. If Lp is the electromagnetic power 
of the pulsar and Lp is the time average (over several months) of this quantity, 
then the pulsed luminosity is L op t = ijLp, where tj is the efficiency, while the SNR 
luminosity is Lsnr = f • Lp + Lo(t), where 0 < / < 1 and the last term (Lo > 0) 
represents the luminosity the remnant would have if the pulsar had no power. At 
maximum brightness of the optical pulses rj > 3-10 _2 /(Xp/Xp)(Afl/47r). The large 
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value of the numerical coefficient constitutes the “optical luminosity problem.” 

Below, we find ri< 10 3 . This implies that emission from the pulsar is beamed 
(AQ << 47 r), or the pulsar wind power is only sporadic ( Lp « Lp), or most 
( Lp — fLp ) of the pulsar spin-down power is either converted into kinetic en- 
ergy of the nebula or reradiated at unobserved frequencies, (or all of the above). 
At any rate, we conclude that the pulsed-beam synchrotron emission model pre- 
sented below can account for all observations if the relatively modest requirement 
f(LpfLp)(ACl/47r)< 10 -1 is met. 

The cavity model is discussed in Section 3, while the constraints implied by 
the data on SN 1987A are considered in Sections 4 and 5. 

2. Difficulties of magnetospheric models 

Optical pulses from the Crab pulsar can originate in that neutron star’s (outer) 
magnetosphere. But if the neutron star in SN 1987A is a weak-magnetic-field (B* < 
10 9 G) “millisecond” rotator (Kristian et.al. 1989, Pacini, Bandiera and Salvati 
1989), it is hard to understand how the optical pulses could arise by an analogous 
process in its magnetosphere. 

Because the Crab pulsar spin rate 27r/Pcrab ~ 200 s —1 ~ 60 times less than 
that of the 1987A neutron star, the emitting area (at the light cylinder radius) can 
be ~ (60) 2 times larger. In additon, the pulsed optical luminosity is an order of 
magnitude smaller in the Crab. The needed Crab optical brightness temperature 
is then ~ 10 6 eV, a value generally exceeded for synchrotron radiation of e* pairs 
created by 7-rays in the outer magnetosphere (Cheng, Ho and Ruderman 1986). 
Such emission mechanisms do not work for the pulsar in SN 1987A for two reasons. 

i) A 10 GeV electron would give peak synchrotron radiation at photon energies above 
100 MeV in the pulsar’s magnetospheTic field. The fraction of energy emitted into 
the optical band would then be very small, ~ 10 -5 of the total radiated synchrotron 
power. 

ii) The detected neutrino burst confirmed that the neutron star in 1987A was formed 
hot, as expected (Hirata et al. 1987, Bionta et al. 1987). The present surface 
temperature of the star should be about 5-10 6 K. The whole magnetosphere between 
the surface of the star and the “light cylinder” (at rj c = P/2n = 3 • 10 6 cm) should 
then be suffused with keV X-rays. In this (black body) X-ray flux, the mean free 
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path for inverse Compton scattering by GeV electrons is ~ 10 3 cm<< r/ c . Therefore 
effective potential drops along the field lines are limited to AU ~ 10 9 V by pair 
plasma created by the Comptonized photons: e + X — » e + 7 followed by 7 -f X — ► 
e + + e - . On the other hand, magnetospheric currents cannot give magnetic fields 
exceeding that of the neutron star. This limits the current flow density along 
open field lines to the Goldreich- Julian value j max — (27r|B|) -1 (ft • B)B (Goldreich 
and Julian 1969), where |Q| = 2n/P. The maximum power of those currents is 
L c = j maz RlQ.~ l AU. Clearly L c > L opt is needed, as the electrons cannot radiate 
more energy than they carry. For L op t = 3 • 10 36 erg/s, a minimum potential drop 
along B of A 17 > 10 14 V is required. This last value is hugely in excess of the 
10 9 V value sustainable without electron pair avalanching. The magnetospheric 
accelerator would thus have been quenched long before it attains the required power. 

It has also been suggested that the neutron star in SN 1987A is vibrating 
with the 0.5 ms period. Wang et al. 1989 proposed cyclotron radiation (in a 
B m ~ 10 12 G magnetic field) of ions powered by surface-penetrating shock waves as 
the mechanism for optical emission. However, it has not been shown how shocked 
ions could gain the necessary velocity perpendicular to B without being fragmented. 
Nor has it been shown how stellar vibration of reasonable amplitude could give rise 
to rapidly recurring shocks of requisite energy. 

We conclude that an origin from within the stellar magnetosphere for the op- 
tical pulsations from SN 1987A has not been plausibly demonstrated for either the 
vibrational or the rotational model. 

3. Pulsar cavities in supernova remnants 

Fax beyond the light cylinder of a pulsar in a vacuum, the spin-down power is carried 
largely in two forms (Rees and Gunn 1974, Kundt and Krotscheck 1977, Kennel and 
Coroniti 1984): 

a) an ultrarelativistic wind, 

b) electromagnetic (EM) fields of the magnetic dipole radiation (from the per- 
pendicular component of the pulsar dipole) and a possible toroidal magnetic field 
(from the spin-aligned part of the dipole) carried with the wind. 

Most of the wind energy is probably due to acceleration of by the very 
strong (time dependent) fields near the pulsar. For a rotating neutron star with a 
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non-spin-aligned dipole moment the pulsar spin frequency would be impressed on 
the electron wind when the electrons are ejected (in a particular direction) from the 
outer magnetosphere and when they are subsequently accelerated. The resulting 
e ± bunch structure would repeat at any (distant) point at the period P of the 
pulsar dipole radiation. If a similar electron injection and wind creation process 
were operative in a strongly pulsating neutron star a modulation at the vibration 
frequency of the magnetic dipole would also be expected. 

When the pulsar is contained within a young SNR the large pressure from the 
pulsar wind and the radiation will create a “cavity” within the remnant. The pulsar 
cavity is terminated by a shock at radius d well within the outer nebula radius D. 
When pulsar emission is the main source of nebular power (Rees and Gunn 1974) 

(d/D) 2 ~ (7 ~ ( 'd/c ) , 

where o is the ratio of the pulsar outflow magnetic energy to the total energy 
density of the wind. For the Crab, Kennel and Coroniti obtain a ~ 3 • 10 -3 and 
^Crab ~ 3-10 17 cm, Kundt and Krotscheck find a ~ 1 and dcrab ~ 10 18 cm. Adopting 
similar values of o for SN 1987A one would then infer a cavity radius d ~ 10 15 cm 
in that SNR, smaller than that in the Crab by roughly the ratio of the SNR ages. 
We do not expect this estimate to be accurate for such a young remnant. However, 
our model only requires that a cavity with radius d < D exist; for SN 1987A, 
D ~ 10 16 cm at the epoch of interest (Papaliolios et al. 1989). 

The outflowing ultrarelativistic bunches of e ± do not radiate significantly in 
the nearly comoving EM waves. To the extent that EM energy is backscattered 
at the cavity wall, they will, however, pass through a magnetic field which may be 
taken to be comparable with that of the preshock incident magnetic field 

B ~ Rem ~ ( aLpjcd 2 ) 1 ^ 2 ~ 2 • 10 -2 (L 39 a- 2 di 2 ) T , 

This value of Rem is similar to the one needed to understand the soft X-ray excess 
emission from SN 1987A, if one assumes equipartition in the nebula (Pacini 1989). 
If u>b = eB /me > 2tt/P, the e ± wind will lose energy in the cavity mostly by 
synchrotron radiation. Had wg < 2ir / P the dominant loss mechanism would have 
been inverse Compton scattering. 
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4. Pulsed emission from the SN 1987A 

In a B ~ 10 -2 G cavity field, the characteristic synchrotron emission frequency is 
~ 10 16 q^Hz, giving optical radiation if 76 = 7/IO 6 ~ 1/6- The fraction of beam 
energy converted to such radiation in a d = 10 15 cm cavity is rj = ju i 2 B {e 2 /mc A )d ~ 
10 — 4 for the same values. Because the optical radiation is emitted almost exactly 
radially, to a distant observer the radiation would appear to be coming from the 
pulsar itself. Thus, cavity and beam parameters of Section 3. could easily give 
the kind of optical luminosity observed from SN 1987A if the wind power were 
~ 10 40 x (Afl/47r) — about ten times the spin-down power of the Crab pulsar 1 if 
emission is isotropic. 

Almost all of the beam power would ultimately be dissipated beyond the cavity 
boundary shock in the surrounding nebula where B is expected to be ~ 10 2 times 
larger than in the cavity. Refer to Section 1. for a discussion of how the current 
upper limit on the bolometric luminosity of the nebula can be satisfied. 

We must now ask what constraints axe imposed on the model parameters by 
insisting that the observed optical (or near infrared) synchrotron light is pulsed 
with the wind frequency 1/P. As shown in the next section, this approach yields 
for the various parameters values close to the ones adopted directly above. We find 
that the size of the nebula places an upper bound 10~ 3 on the efficiency of 

radiation allowed by the model. 

A critical assumption is that the relativistic electrons synchrotron radiate in 
an ordered EM field of wavelength cP. This guarantees that the deflection from 
the radial direction of the radiating never exceeds am angle (#0, eq. [P 10]) less 

1 The expected pulsed cavity emission from the Crab can be scaled from that 
from SN 1987A. For the “optical” frequency oJCrab /^1987 = [7 2 i?]crab/[7 2 #]i987 ~ 
[^VoT7p/d\cr&b/[^V^Lp/<I\\9»i- For comparable 7 and oLp, u>c r ab ~ ^1987/500 
or A(Crab) ~ 10 2 /xm. With similar approximations and assumptions the ratio of 
pulsed cavity emission luminosities from the Crab and SN 1987A is the ratio of 
the values of crLpj/d, again correxponding to a reduction of about 500. Thus, 
the Crab’s pulsed cavity far IR luminosity would be ~ 10 33 erg/s. A bump of 
about this magnitude appears in the near (A < 3.5/jm) IR pulse shape of the Crab 
(Middleditch, Pennvpacker and Burns 1983). 
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than the critical one beyond which the pulses would be washed out. If, instead, 
the field had been a collection of randomly oriented domains of size cP the average 
total deflection would have been too large, 80 (d/cP) 1 / 2 . 

5. Constraints on pulsed beamed synchrotron emission implied by the 
SN 1987A data. 

By assumption, the optical signal is due to synchrotron radiation of relativistic e ± 
(energy 7 me 2 ) in transverse magnetic field of strength B alternating in direction 
with wavelength cP. Before entering an assumed emission zone of radial extent 
1 , the electrons travel radially outwards a distance d — l from the neutron star. 
The electrons radiate into a narrow forward cone of apex angle ~ lfy about their 
instantaneous velocity direction, which is, itself, at an angle to the initial (radial) 
direction of flight. The latter angle is not greater than some maximum deflection 
angle 6q (eq. [P 10 ]). Thus, the cross-sectional area of the emission region seen by 
an observer at infinity is fs 7 rt 2 , where 

b « d8, (PI) 

and 

e « 8 0 + 1/7 « 1 . (P2) 

For the purposes of computation we take the optical brightness temperature to 
be kTb — 10 3 GeVx (fc 2 / 10 12 cm 2 ) 1 and the synchrotron frequency to be 

y 2 eB/mc = 2 eV/h, (P3) 

(i.e. t 2 P/ 10 8 G = 2 ) to obtain (nearly) optimum efficiency of optical detection. 
The maximum extent of the nebula, D ~ 10 16 cm places an upper bound on the size 
of the emitting region and its radial distance from the star: l < D, d < D. 

Notation 

7 — initial Lorentz factor of the radiating electron 

8 — maximum angle between line of sight and initial direction of electron motion 
8 0 — maximum deflection angle of electrons 

d — radial distance from the neutron star to the emitting region 
/—radial extent of the emitting region 



■Kb 2 — area of emission seen by observer 

At — maximum allowed differential time of arrival (t.o.a.) 

a, /?, A, di 6 < 1 — dimensionless parameters not greater than unity 

r, G > 1 — dimensionless parameters greater than unity 

THE CONSTRAINTS 

A class of constraints is introduced by the requirement that the optical pulses 
not be washed out. Let the upper bound on the differential spread in time of arrival 
of all photons in a pulse be At = aP / 5 = a x 10 -4 s, i.e., 

cAt = a x 10 6 ' 5 cm, a < 1, (P 4 ) 

Any initial spread in energies (me 2 A7) of e ± leads to a constraint l^< ^(cAt) /A7, 
less stringent than the following. We define 

C=i(* , r I + l)«l+i r ! Dj + * (P5) 

The differential t.o.a. constraint from time of flight delay of the emitting e ± gives 

1 = XG~W(cAt), A < 1, (P6) 

Differential t.o.a. because of different path lengths due to the transverse extent of 
the emitting region gives b = /38~ 1 (cAt), and hence 

d = /36~ 2 (cAt), 0<1, (P7) 

(Strictly speaking A -j- (} < 1, but we are not concerned with factors of 2.) We note 
the following limits: 

8 0 « 1/7 => G — 1, 

80 ~ 1 / 7 = 4 >( j ~ 2 , 

< 9 0 » 1/7 =►<?= 

The inferred brightness temperature places a lower boimd on the electron energy 

7 = 10 5 - 3 r y9“ 2 a- 2 0 2 , r > 1. (PS) 
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The efficiency of conversion of the electron energy to optical is rj ~( Synchrotron 
power)x(ymc 2 )- J X l/c, i.e., 

1 - §7 X 10“ («) 

where eq. (P3) was used to eliminate B. 

Since the magnetic field traversed by the e* alternates in direction, the appro- 
priate expression for the deflection angle is #o ~ PeB /( 27 rjmc), i.e. 
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11.4 


(P 10 ) 


while the size of the nebula limits d and l, 

l < d = d 16 x 10 16 cm, d 16 < 1. 


(P 11 ) 


RESTRICTIONS ON 7 

Consider the constraints (P6) to (P9) in the following two cases. 

0#o si-' 

From eq. (P10) this regime holds iff 7 2 ^ 10 11 * 4 , i.e. 7 >, 10 5-7 . However, the 
constraint (P6) with the subsidiary condition (Pll) then gives a low efficiency, since 
l = aA-yg x 10 18 ' 5 cm. Hence, orA < 10 _2 <fi6 and therefore r ]< 10 -4 di6. 

ii)#o >> 7 -1 

From point i) above, this can only hold if 7 < 10 5 ' 7 . Now, from eqs. (P5) and 
(P10) 

«»( 2C?)V7~ 10 114 / 7 3 . 

Hence, using also eq. (P8), 7 = x 10 4 ° < 10 5 7 , i.e., 

10 4 < 7 < lO 5 ' 7 . (P12) 

Now, from eqs. (P6), (P7) and (P9), respectively, 

1 « (oA) 74 6 x 10 8 '°cm « (A//?)r*(<*£)- 5 / 7 x 10 8O cm, (P13) 

d « (aft) 7® x 10 ' ''cm ^ r?(ar/?) - 7 x 10 7 ' 7 cm, 
r? * (oA) 74 3 x 10~ 6 4 « j f3)T 2 ^ [a j3)^ x 10“ 6 - 4 . 
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(P14) 

(P15) 


Recall that a is the precision with which fine structure can be observed in the optical 
pulses in units of P/5 = 0.1 ms and di 6 is the maximum distance of the emitting 
region from the pulsar in units of 10 16 cm. In principle, a is an observable number. 
The free parameters A, f3 (and T) were introduced to replace upper (lower) bounds 
with equalities. 

To satisfy l < d we must have A < 0/2 < 10~ 0 3 . But (P13) and (Pll) require 
(aA)7| < 10 8 0 di6, giving, upon substitution in eq. (P15), rj < 10'' 2 ' 4 (aAdi 6 ) 1 / 2 . 

We conclude that this model allows a maximum efficiency of 

W^S-lO-Vae) 1 ' 2 , (P16) 

occurring for 7 = 10 5 ' 5 (die) 1 / 6 and the most favorable values possible of the free 
parameters (/? = 1, A = 0.5). 
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